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ABSTRACT 

Knowledge of mass and concentration of galaxy clusters is crucial to understand their forma- 
tion and evolution. Unbiased estimates require the understanding of the shape and orientation 
of the halo as well as its equilibrium status. We propose a novel method to determine the 
intrinsic properties of galaxy clusters from a multi-wavelength data set spanning from X-ray 
spectroscopic and photometric data to gravitational lensing to the Sunyaev-Zel'dovich effect 
(SZe). The method relies on two quite non informative geometrical assumptions: the distribu- 
tions of total matter or gas are approximately ellipsoidal and co-aligned; they have different, 
constant axial ratios but share the same degree of triaxiality. Weak and strong lensing probe 
the features of the total mass distribution in the plane of the sky. X-ray data measure size 
and orientation of the gas in the plane of the sky. Comparison with the SZ amplitude fixes 
the elongation of the gas along the line of sight. These constraints are deprojected thanks 
to Bayesian inference. The mass distribution is described as a Navarro-Frenk-White halo 
with arbitrary orientation, gas density and temperature are modelled with parametric profiles. 
We applied the method to Abell 1689. Independently of the priors, the cluster is massive, 
M20Q — (1-3 ± 0.2) x 1O 15 M , and over-concentrated, C200 = 8 ± 1, but still consistent 
with theoretical predictions. The total matter is triaxial (minor to major axis ratio ~ 0.5 ±0.1 
exploiting priors from TV -body simulations) with the major axis nearly orientated along the 
line of sight. The gas is rounder (minor to major axis ratio ~ 0.6 ± 0.1 ) and deviates from 
hydrostatic equilibrium. The contribution of non-thermal pressure is ^20-50 per cent in inner 
regions, < 300 kpc, and ~ 25 ± 5 per cent at ~ 1.5 Mpc. This picture of A1689 was obtained 
with a small number of assumptions and in a single framework suitable to application to a 
large variety of clusters. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: Abell 1689 - galaxies: 
clusters: intracluster medium - methods: statistical 



1 INTRODUCTION 

Clusters of galaxies, the most recent bound structures to form in 
the Universe, are excellent laboratories for precision astronomy 
jVoIt||2005) . Their use in cosmological tests relies on accurate 
measurements of their mass and concentration (Meneghett i et al.| 
|2010| [Rasia et a l.|2012| l. Assessing the equilibrium status is also 
crucial in determining evolution and mechanisms of interaction 
of baryons and dark matter ( |Lee & Si ito 2003 , Kazantzidi s et al.| 
|2004[ >. An unbiased look at the cluster properties must take into 
account their shape and orientation too jOguri et al.||2005") . The 
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intrinsic form shows how material aggregates from large-scale per- 
turbations ( |West|1994||Jing & Suto|2002| >. On the other hand, esti- 
mations of mass, inner matter density slope and concentration may 
be biased if derived under the assumption of spherical symmetry 
( |Gavazzi|20"03||Meneghetti et al.|2010{|Rasia et al.|2012) . 

Assessing the intrinsic shape and orientation is also criti- 
cal when comparing observations with theoretical predictions. The 
slope of the concentration-mass c(M) relation of galaxy clusters 
is consistent with TV-body simulations, though the normalization 
factor is higher jComerford & Natarajan|2007||Ettori et al.|20i0| >. 
Galaxy clusters selected according to their gravitational lensing 
strength or X-ray flux may form biased samples (Meneghetti et al. 
|201 1| >. Triaxial halos are more efficient lenses than their more 
spherical counterparts ( Oguri & Blandford 2009 ) with the strongest 
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lenses in the Universe expected to be a highly biased population 
preferentially orientated along the line of sight. In fact, the over- 
concentration problem is less prominent in analyses accounting for 
triaxiality jOguri et al.|2005[|Sereno & Zitrin|2012| >. 

We can test the intrinsic three dimensional form of a popu- 



lation of astronomical objects with statistical approaches l Hubble 



1926, Noerdlinger 1979, Binggeli 1980, Binney & de Vaucouleurs 
19811 |Fasano & Vio||1991| |de Theije et al.||1995| |Mohr et ~ 
1995||Basilakos et al|2000| |Cooray|2000| |Thakur & Chakraborty 



200Tj|Alam & Ryden|2002] Ryden|1996||Plionis et al.|2004||Paz 



et aH2 006; Kawahara |2010) . Clusters of galaxy can be observed 
with very heterogeneous data-sets at very different wave-lengths 
from X-ray surface brightness and spectral observations of the 
intra-cluster medium (ICM), to gravitational lensing (GL) observa- 
tions of the total mass distribution to the Sunyaev-Zel'dovich effect 
(SZe) in the radio-band. This enables us to tackle the structure of 
a cluster with a multi-probe approach (|Zaroubi et al.| 1998 Reblin-| 
|sky|2000||Dore et al.|200Tj|Puchwein & Bartelmann|2006| l. Only 
a few works have tried to infer shape or orientation of single ob- 
jects. Combined use of X-ray and SZe data enables to constrain 
the shape of the ICM without any assumption regarding equilib- 
rium. jDFTjjippjFerar] ([2005} and |Sereno et al.| ( |2006) first stud- 
ied a sample of 25 clusters finding signs of a quite general triaxial 
morphology. Mahdavi & Chang (201 1 1 constrained the minimum 
line-of-sight extent of the hot plasma of the Bullet cluster with a 
model-independent technique. 

Lensing observations offer an alternative way to study triaxial- 
ity. Projected mass distributions from either weak or strong lensing 
can be deprojected exploiting some a priori assumptions on the in- 
trinsic ^shapejiJOgumjnir]2005^ |Morandi et al.| 



|20TT| [Sereno et al.|2010||010b||Sereno & Umetsu|2011 

Abell 1689 (A1689) is a very luminous cluster at redshift 
z = 0.183 ( |Broadhurst eTaI1[2005l |Limousin et al.||2007> . The 
mass distribution in the inner < 300 kpc regions of the cluster 
has been accurately determined by strong lensing analyses that 
favoured a quite concentrated mass distribution (Broadhurst et al. 
|2lX)5l|HaIkbra et al.|2006"l|Limousin et al.|2007||Coe et al.|2010| >. 
On the larger virial scale, different weak lensing analyses suggest 
somewhat different degrees o f concentration ^Umetsu & Broad 



hurst|2 008 , Ume tsu et al."|200"9"l|Corless et al.|2009HLimousin etal. 



2007| >. A further puzzle is the conflict between X-ray and lensing 
analyses, with lensing masses exceeding estimates derived under 
the hypothesis of hydrostatic equilibrium by 30-40 per cent in the 
inner regions ( |Peng et al.|2009| l. 

A1689 has been object of a number of triaxial analyses jOguri| 
|et al.|2005||Corless et al.|2009) . |Sereno & UmetsulpOTT) analyzed 
weak and strong lensing data and found evidence for a mildly tri- 
axial lens (minor to major axis ratio ~ 0.5 ± 0.2) with the ma- 
jor axis nearly aligned with the line of sight. The halo was over- 
concentrated but still consistent with theoretical predictions. |Peng| 
et al. (20091 recognized that a prolate configuration, aligned with 
the line of sight and with an axis ratio of ~ 0.6, could solve the cen- 
tral mass discrepancy between lensing and X-ray mass estimates. 
|Morandi et al.| l |2011} combined lensing and X-ray data assuming 
the cluster to be aligned with the line of sight and fixing the rela- 
tion between gas and matter shape. They found an axial ratio for 
the matter distribution of ~ 0.5. |Sereno et al.| ( [2012fr combined X- 
ray and SZe data to model the gas shape and orientation by using 
Bayesian inference. Their analysis favoured a mildly triaxial gas 
distribution with a minor to major axis ratio of 0.70 ± 0.15, pref- 
erentially elongated along the line of sight, as expected for massive 
lensing clusters. 



Here, we propose a method to infer the intrinsic shape and 
orientation of clusters based on weak and strong lensing (WL and 
SL) observations plus deep X-ray and SZe observations. Lensing 
gives a picture of the total projected mass distribution without any 
assumption on the equilibrium status of the cluster. X-ray plus SZe 
observations fix the size of the gas distribution along the line of 
sight and in the plane of the sky ( Sereno 2007 , S ereno et al.|2012| . 
The method exploits Bayesian inference to study a number of vari- 
ables larger than the number of observational constraints. This en- 
ables us to investigate both intrinsic shape and orientation. The ver- 
sion of the method detailed in the present paper joins and integrate 
the lensing analysis in|Sereno & U metsu (201 1 ) with the X-ray plus 
SZe investigation of |Sereno et al. 1 2012) . 

The paper is organised as follows. In Sec. [2] we discuss how 
triaxial ellipsoids project and the relation between total matter and 
gas distributions. Sec. [3] summarizes the X-ray plus SZe analysis. 
Secsj4]and|5]are devoted to the WL and SL parts, respectively. In 
Sec. [6] we show how to combine the different data-sets to infer the 
intrinsic properties. Results are discussed in Sec. [7] The hydrody- 
namical status of the cluster is wrote about in SecflS] Comparisons 
with previous works are listed in Sec. [9] Final considerations are 
contained in Sec. 1101 

Throughout the paper, we assume a flat ACDM cosmology 
with density parameters Qm =0.3, Qa = 0.7 and Hubble constant 
100/1 kms _1 Mpc 



Hi 



Ii 



0.7 ± 0.014 (Komatsu et al 



2011 1. At the A1689 distance, 1" corresponds to 2.15 kpc//i (= 



3.08 kpc). 



2 TRIAXIAL HALOS 



In this section we describe how we model the total matter and the 
gas distribution and how we relate them. 



2.1 Matter profile 

The ellipsoidal Navarro-Frenk-White (NFW ) density profile, 

Ps 



Pnfw = 



(C/V s )(l + C/r s ) 2: 



(1) 



provides a very good fit to the density distribution of dark matter 
halos in high resolution TV-body simulations (Navarro et a l.|1996| 
[7997l|Jing & Suto|2002} . In Eq. (]}, ( is the ellipsoidal radius. The 
shape of the halo is determined by the axial ratios, which we denote 
as qi (minor to major axial ratio) and q2 (intermediate to major 
axial ratio). The eccentricity is et — ^Jl — qf. The orientation 
of the halo is fixed by three Euler's angles, 8, ip and ip, with i9 
quantifying the inclination of the major axis with respect to the line 
of sight. 

The NFW density profile can be described by two parame- 
ters, the concentration and the mass. By definition, r2oo is such 
that the mean density contained within an ellipsoid of semi-major 
axis r2oo is 200 times the critical density at the halo redshift, 
Per <|Corless et al.||20"0"9l |Sereno et aT"||2010[ |Sereno & Umetsu| 
2011 1. Then, the corresponding concentration is C200 = >"20o/Vs. 



M200 is the mass within the ellipsoid of semi-major axis r2oo, 
M200 = (8007r/3)(?ig2rlooPcr. 

The projection into the sky of the ellipsoidal 3D NFW halo 
is an elliptical 2D profile l |Stark|1977l|Sereno|200"7j|Sereno et al.| 
|010b[ >. The convergence k, i.e., the surface mass density in units 
of the critical density for lensing, E cr = (c 2 D B )/(4irG Dd Da s ), 
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where D s , Dd and Dds are the source, the lens and the lens-source 
angular diameter distances respectively, can be written as 



knfw(s) 



1 



= arccosh I - 



x is the dimensionless elliptical radius, 

x = ar sP , Z=[x 2 +x 2 2 /(l~e) 2 )} 1/2 , 



(2) 



(3) 



where e is the ellipticity and xi and X2 are the abscissa and the 
ordinate in the plane of the sky oriented along the the ellipse axes, 
respectively. 

The strength k s and the projected length scale r s p are related 
to mass and concentration and depend on the shape and orienta- 
tion parameters too. Explicit formulae can be found in Sereno et al. 
pOlfj) . 

Since we can only observe projected maps, the problem of 
determining the 3D orientation and shape of halos is intrinsically 
degenerate. Ellipsoids project into ellipses (Stark 1977 1. Even with 
an ideal multi-probe data-set without noise, we can only measure 
three observable quantities which help us to constrain the five un- 
known intrinsic properties (two axial ratios and three orientation 
angles) of the ellipsoidal halo (Sereno 20071. The three measurable 
quantities are the ellipticity e, the orientation 9 t and the elongation 

The parameters e and 9 t characterise the projected ellipse in 
the plane of the sky. The ellipticity e gives a measurement of the 
width in the plane of the sky. It is defined as 1 — b p /a p , where 
bp and a p are the minor and major axis of the projected ellipse. 6 ( 
measures the orientation in the plane of the sky of this ellipse. 

eA quantifies the extent of the cluster along the line of sight. It 
is the ratio between a p and the size of the ellipsoid along the line of 
sight ijSereno et al. 2012. see figure 1). The smaller eA, the larger 
the elongation along the line of sight. If eA < 1, the cluster is more 
elongated along the line of sight than wide in the plane of the sky. 

These three observable quantities depend on the intrinsic axial 
ratios and on the Euler's angle (Binggeli 1980, Sereno 2007). 

2.2 ICM 



Observations (jKawahara 2010) and theory JBuote & Humphrey| 
|2012| l suggest that also the density of the intra-cluster medium 
(ICM) is nearly constant on a family of similar, concentric, coax- 
ial ellipsoids. Modelling both the gas and the matter distribution as 
ellipsoids with constant eccentricity is formally wrong in haloes in 
hydrostatic equilibrium. If the cluster of galaxies is in hydrostatic 
equilibrium the gas distribution traces the gravitational potential. 
Given an ellipsoidal gas density, the gravitational potential is el- 
lipsoidal too and can turn unphysical for extreme axial ratios, with 
negative density regions or unlikely configurations. On the other 
hand, even in hydrostatic equilibrium, the ellipsoidal approxima- 
tion for the gas is suitable in the inner regions or when small ec- 
centricities are considered. 

If the potential is ellipsoidal, the matter distribution that origi- 
nates it can not be ellipsoidal. However, dark matter haloes formed 
in cosmological simulations typically have radially varying shapes 
too {Kazantzidis et al. 2004 1, which might be compatible with an 
ellipsoidal potential. 

If the matter halo isodensity surfaces are triaxial ellipsoids, the 
isodensity surfaces of the intracluster gas are well approximated as 
triaxial ellipsoids with eccentricities slowly varying with the ra- 
dius ( |Lee & Suto 2003 2004). The ratio of eccentricities of gas 



with e5 CM /ef at ~ 0.7 for i = 1, 2. Furthermore, the variation in 
eccentricity is usually smaller than the observational error on the 
measured ellipticity of the X-ray surface brightness map. 

A further complication is that microphysical processes such as 
radiative cooling, turbulence or feedback mechanisms strongly af- 
fect the shape of the baryonic component and make the ICM shape 
more triaxial and with a distinctly oblate shape towards the cen- 
tral cluster regions compared to the underlying dark matter poten- 
tial shape jLau et al.|201 1| >. The gas traces the underlying gravita- 
tional potential better outside the core, even if radiative processes 
can make the ICM shape rounder. These mechanisms can be effec- 
tive so that assuming that the overall triaxiality of the gas is strictly 
due to the underlying shape of the dark matter potential can be mis- 
leading. 

Finally, density and shape of the gas distribution have a strong 
correlation, whereas temperature is essentially uncorrelated l [Sams-| 
|ing et al.|2012) . Then the proper modelling of a varying eccentricity 
would require an independent measurements of the density profile. 

Under these circumstances, we can make two simplifying but 
quite non-informative working hypotheses to relate gas and total 
matter distributions. Firstly, we conservatively model both the to- 
tal matter density and the gas distribution as coaligned ellipsoids 
with fixed, but different, axial ratios. Secondly, we set the two dis- 
tributions to have the same triaxiality parameter T(= e 2 /e 2 ).lf 
two distributions have the same triaxialities, then the misalignment 
angle between the orientations in the plane of the sky is zero ( Ro- 
manowsky & Kochanek 1998 1. This is in agreement with what ob- 
served in A1689, where the centroid and the orientation of the sur- 
face brightness map l |Sereno et al.||2012) coincide with those of 
the projected mass distribution as inferred from lensing ( |Sereno &| 
|Umetsu|20iT) . 

The axial ratios of the gas distribution, q} CM , can be expressed 
in terms of the corresponding axial ratios of the matter distributions 

as 



^Jl-(e lCM /e M ^) 2 {l-qf). 



(4) 



= T ICM , then ei CM /er at 



ICM /„Mat 



7 



Being 7~ Ma 

g ICM ^gMat 

ial ratios to three: q± and q 2 for the matter and g} CM for the gas 
q\ cm is determined by T(qi , q 2 ) and q\ CM , 



The above assumptions limit the number of free ax- 



ICM i 

q 2 =4/1 



1 - b\ 



ICMJ2 



(5) 



(e ) and matter (e a ) is nearly constant up to the length scale, 



In a triaxial analysis, results strongly depend on the minor to ma- 
jor axial ratio whereas the inference about the intermediate axis is 
more affected by priors. Then, using q\, q 2 and q\ CM as free pa- 
rameters is not really a limitation with respect to considering all 
four axial ratios as free. Finally, as expected from both theory and 
observations, we required that the gas distribution is rounder than 
the matter one, q?" ^q{ CM . 

Under these hypotheses, total matter and gas have different el- 
lipticities and elongations but share the same orientation 8 € . Since 
we model the gas isodensities as ellipsoids, they project into el- 
lipses in the plane of the sky. The same degeneracy plaguing the 
inference of shape and orientation of the matter distribution affects 
the gas too. 
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Figure 1. Probability distribution of the elongation of the gas distribution 
as inferred from the X-ray plus SZe analysis. 



To asses realistic probability distributions for the parameters 
we performed a statistical Bayesian analysis. The Bayes theorem 
states that 

p(P|d) cx £(P|d)p(P), (8) 

where p(P|d) is the posterior probability of the parameters P 
given the data d, £(Pjd) is the likelihood of the data given the 
model parameters and p{P) is the prior probability distribution for 
the model parameters. 

The posterior probability distribution for the elongation of the 
gas eJf M is plotted in Fig. [I] The elongation together with the pro- 
jected ellipticity and orientation enables us to constrain the intrinsic 
shape of the gas. Whereas the projected ellipticity and the orienta- 
tion angle of the projected isophotes of the gas can be precisely 
determined with X-ray data alone, the measurement of the elonga- 
tion requires a combined X-ray plus SZe analysis. 



3 X-RAY AND SZE 

The shape and orientation of the gas distribution can be inferred 
with a combined analysis of X-ray and SZe data (Seren o~et al.| 
|2012| >. The first set of constraints comes from the analysis of the 
projected maps. The ellipticity and the orientation of the isophotes, 
which could be measured with either X-ray or SZe alone, are re- 
lated to the size of the gas distribution in the plane of the sky. The 
inference of the elongation along the line of sight requires a com- 
bined analysis. This can be done thanks to the different dependence 
of X-ray and SZe observables on the gas density. 

The analysis of X-ray and SZe data of A1689 was performed 
in |Sereno et al.| ([2012l. Here, we summarize the main results and 
what is needed in the present paper. We refer to |Sereno et al.|p012) 
for details. Images in the 0.7-2 keV band taken by Chandra are 
well approximated by a series of concentric, coaligned ellipses with 
e x = 0. 15±0. 03 and orientation angle 9? = 12±3deg (measured 
North over East). The 3D electronic density and temperature were 
modelled with the following parametric profiles ( Vikhlinin et al.| 
|2006l|Ettorietal.|2009| >, 



n c = n 



1 + 



-3/3/2 



1 + 



(6) 



where no is the central electron density, r c is the core radius, r t (> 
r c ) is the tidal radius, /3 is the slope in the intermediate regions, and 
7 is the outer slope. For the temperature profile, we used 



T = 



To 

i + (CAt) 2 ] " 



(7) 



where To is the central temperature and the radius n describes a 
decrement at large radii. The parametric forms in Eqs. |6]|7]l were 
motivated by the absence of a cool core. 

Profiles based on Eqs. l(6]|2) were then compared to observa- 
tions of X-ray surface brightness (SB), spectroscopic temperature 
and SZ decrement. The surface brightness observed by Chandra 
was collected in 68 elliptical annuli up to £ < 1200 kpc. The tem- 
perature profile from the XMM satellite was binned in 5 elliptical 
annuli up to £ < 900 kpc. We considered the integrated Comp- 
ton parameters Y within the circle of radius r2soo ~ 600 kpc 
measured from 7 observatories (BIMA, OVRO, AMIBA, SuZIE, 
WMAPS, SZA, SCUBA). We re-did the analysis as in |Sereno et al~| 
l |2012) , but we used an improved version of the code, with some 
better numerical algorithms for numerical integration. 



4 WEAK LENSING 

We fitted the surface mass density of A1689 with a projected ellip- 
soidal NFW halo. The weak lensing convergence map, described 

, was obtained from 



as "2D MEM" in |Umetsu & B roadhurst 
Subaru data and covers a field of ~ 30' x 24' (21 x 17 grid pix- 
els with pixel size of 1.4') ( jUmetsu & Broadhurst|2008l 



Umetsu 



|et al.|2009] >. We closely followed the method employed in |Sereno 
& Umetsu ( ]201 1| >, which we refer to fo r more details. The likeli - 
hood can be written in terms of a Xwl 1 Oguri et al.|2005 2010 1, 
with 



2 

Xwl 



[ K °bs(r0 - re(ri)] (V *) [Kobsfo) - wfa)] (9) 



where K, b s is the measured convergence map and V -1 is the in- 
verse of the pixel-pixel covariance matrix. 

One major difference with Sereno & Umetsu (201 1 1 is in the 
priors used to determine the NFW projected parameters. Whereas 
in |Sereno & Umetsu] ( |201 1| > uniform priors were employed, here 
we exploit constraints from the X-ray analysis on centroid position 
and ICM orientation. The position of the centroid of the X-ray sur- 
face brightness, {a, 8} = {197.87274, -1.3400533} deg with an 
accuracy of ~ 1.2", is consistent with the centre of the BCG and 
with the centroid of the total matter distribution as estimated from 
lensing (Sereno & Umetsu 201 1 1. The orientation of the X-ray map 
is consistent too with the orientation of the matter halo as estimated 
from lensing. The posterior probability is then 



oc exp(- XW L/2) x Af(9 e ;9f, 59?) 
x Af(0i, o ; 9?, 69?) x Af(9 2 , ; 9%, 89%) 



(10) 



where N{P; n, o) denotes a normal distribution for the parameter 
P centred in /i and with variance a. 9i,o and #2,0 are the Cartesian 
coordinates of the projected centre of the matter distribution in a 
reference system centred in the BCG galaxy. 9? and 6 2 are the 
observed coordinates of the X-ray centroid. 

The parameter space was explored with Markov chains. Chain 
convergence was checked by requiring that the standard var(chain 
mean)/mean(chain var) indicator was less than 1.2. Results are 
summarised in TableQ] 

Analyses based on WL alone are not very effective in deter- 
mining the projected ellipticity. The relative constraints on the in- 
trinsic axial ratios are not sharp and are dominated by the a priori 
hypothesis employed. Here, since we exploited the X-ray informa- 
tion on the halo orientation in the plane of the sky, the ellipticity is 



© 0000 RAS, MNRAS 000, 000-000 



Multi-probe of A 1689 5 



NFW k s 


W J.,U "z^J 

["] ["] 




0, 
[deg] 


[kpc] 


0.62±0.18 


Weak Lensing 
-0.8±1.2 -3.1 ±1.1 


0.20±0.10 


11.6±2.9 


220 ±50 


0.27±0.02 


Strong Lensing 
-3.0±0.7 -4.8 ±2.7 


0.16±0.03 


13.6±5.6 


910 ±120 


0.35 ±0.01 


Weak plus Strong Lensing 
-3.8±0.5 -0.1±1.0 0.11 ±0.02 


12.7±2.4 


410 ±20 



Table 1. Projected NFW parameters inferred with either the weak or the strong or the combined lensing analysis. The final WL distributions for the centroid 
position and the orientation angle mirror the X-ray derived priors. The orientation angle 6 e is measured north over east and r s p is the projected length scale. 
Convergence is normalised to a reference source redshift (z B = 2.0). Central location and dispersion of the marginalized posterior density functions (PDFs) 
are computed as mean and standard deviation. 



better constrained. With respect to the WL only analysis in Sereno 
|& Umetsu] ( |201 1[ >, the final distributions on the projected parame- 
ters are more peaked and with reduced tails. This is clearly seen in 
the probability function for e, which would be nearly flat without 
the X-ray priors. 



5 STRONG LENSING 

We performed a strong lensing analysis of the inner regions of 
A1689 and obtained a pixelated map of the surface mass density 
by use of the PixeLens software (Saha & Williams 20041. We con- 
sidered multiple image systems with confirmed spectroscopic red- 
shift, whose detailed description can be found in Limousi n et al.| 
( |2007} ; |Coe et al.|(2010[ >. In a second step, we fitted the map with a 
projected NFW profile. The method is described in detail in |Sereno] 
|& Zitrin||20T2) 

PixeLens cannot handle all the multiple image systems of 
A 1689 at once. We then divided the SL systems in two groups and 
analyzed each group separately. The first group includes the sys- 
tems 1, 2, 7, 10, 11, 18, 22, 24, 40 according to the notation in 
|Limousin et al.|p007| >, for a total of 31 multiple images. The sec- 
ond group includes the systems 4, 5, 15, 17, 19, 29, 33 and 35 (28 
multiple images). For each group, we computed 500 convergence 
maps within a radius of 94"around the BCG on a grid of 861 pixels 
with a pixel size of 12.6 kpc/ft. 

As explained in Sereno & Zitrin (2012), we excised a central 
region of 40 kpc//i to minimize the effects of miscentering and 
baryonic physics ( [Umetsu et al.|20lTj . We also excluded the outer 
pixels where the logarithmic density slope artificially takes values 
smaller than -3. 

We then modelled the convergence with a projected NFW pro- 
file. For each convergence map, we looked for the minimum of 



obs(l"i) — KNFw(l"i)] 



(ID 



where the sum runs over the pixels. From the derived ensemble of 
maximum likelihood parameters, we obtained the posteriori distri- 
bution of projected NFW parameters. We repeated the procedure 
for the two groups of images, ending up with results statistically 
not distinguishable. 

Results were consistent with previous SL analyses based on 
the same image systems ( |Limousin et al.|[2007| |Coe et aL]|2010| 
|Sereno & Umetsu 20lT] and references therein). Mass density 
profiles obtained with PixeLens are usually shallower than recon- 



structions based on parametric models JGrillo et al.|201(M Umetsu 
let al.||2012l. In comparison with the mass profiles in Limousin 
|et al.| (|2007 1, who used Lenstool, or |Sereno & UmetsuH2011) , who 
used Gravlens, PixeLens retrieves a profile shallower in the inner 
< 90 kpc//i. On the other hand, the projected masses within the 
Einstein radius (< 300 kpc) are fully consistent. The excision of 
the inner 40 kpc/h in our analysis makes the impact of the differ- 
ences between different methods nearly negligible. 

Results from SL are marginally consistent with the WL anal- 
ysis, see Table [T] WL favors more concentrated profiles but the 
strong degeneracy between the lensing strength k s and the pro- 
jected length-scale r B p makes the two analyses compatible. 

The measurement of the orientation angle with the strong lens- 
ing analysis of the core region is in very good agreement with the 
estimate from the X-ray analysis and further supports the working 
hypothesis that the gas and the total matter distribution share the 
same triaxiality degree and are then co-aligned in projection. 



6 DEPROJECTION 

Results from either gravitational lensing or X-ray plus SZe mea- 
surements can be combined to deproject the observations and infer 
the three dimensional structure and orientation of the matter and 
gas distributions. The combined X-ray plus SZ analysis enabled us 
to infer the width of the ICM in the plane of the sky (parameterized 
in terms of the ellipticity e ICM ) and its size along the line of sight 
(expressed as the elongation e^f M ). 

The lensing analysis describes how the total matter density 
projects in the plane of the sky. The orientation and the ellipticity of 
the projected surface density are related to the intrinsic shape and 
orientation of the total matter halo. How the convergence varies 
with the radius, i.e., the information contained in the parameters 
k s and r s p of the NFW halo, constrains the functional form of the 
density, i.e, the mass and the concentration. 

The orientation and the shape of the matter halo significantly 
affect lensing. The more the cluster is elongated along the line of 
sight, the more the apparent convergence is boosted and the smaller 
the projected length scale in the plane of the sky. 

Information from X-ray, SZe and lensing can be brought to- 
gether. The likelihood function combining the results from the pre- 
vious sections can be written as 



c = r GL x £ ICM 



where £ GL = £ WL x £ SL and C x 



(12) 

£(K s ,r sP ,£,0 e ) 
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l |Sereno & Umetsu|201 1| >. For the likelihood of the X-ray plus SZe 
part i Sereno et al.|2012 1, 



,ICM 



(2 7 r) 1 /2 (Te 
x P(ef M 



exp 



AeT) 



■ exp 



2ol 



1 /Ae^ 2 



(13) 



where P(e£ M ) is the marginalized posterior probability distribu- 
tion for the elongation parameter of the gas obtained in Sec. [3] The 
parameter Ae^ s quantifies the additional unknown statistical and 
systematic uncertainty on the elongation. We let it follow a normal 
distribution (D'Agostini 2003), centred on and with dispersion 



8e s ^ s ~ 0.13 (Sereno et al. 



2012) 



The overall ellipsoidal cluster model describing both ICM and 
dark matter has 9 free parameters. For the total matter halo: 2 pa- 
rameters determining the profile, i.e, the mass A/200 and the con- 
centration C200; 2 axial ratios (qi and 92) and 3 orientation angles 
if and 

For the gas distribution: the axial ratio q\ CM , determining the 
shape of the ICM; 1 parameter, Ae^ s , quantifying the systematic 
uncertainty on the elongation of the gas distribution. Mass and gas 
share the same orientation angles whereas the second axial ratio of 
the ICM is fixed by the remaining parameters thanks to the priors. 

The distributions of matter and gas project into ellipses into 
the plane of the sky. Lensing constraints reduce to four param- 
eters. The analysis of the ICM adds two more measured quanti- 
ties. The intrinsic parameters have to be then determined from a 
smaller number of observational constraints. The problem is under- 
constrained ( Sereno 2007 1. 

From the lensing analysis, we can measure the 4 projected 
parameters of the NFW halo (the lensing normalization re s , the 
projected length-scale r s p, the ellipticity e and the orientation 6 e ). 
From the analysis of the ICM, we can estimate the size of the gas 
in the plane of the sky (e ICM ) and along the line of sight (e^f M ). 
The information on the orientation of the X-ray isophotes in the 
plane of the sky, 6\ CM , was exploited as a prior in the WL analysis 
to better constrain 6 t . The estimate for the systematic error Se s ^ s 
provides another constraint. 

Let us summarize the dependences of each observed quantity. 
The ellipticity e ICM and the elongation e l £ of the gas distribu- 
tions are functions of gl CM , q}p M , and of the orientation angles 
■& and if. Due to the imposed equality of triaxiality between the 
mass and gas distributions, we can express q\ CM in terms of qi CM , 
qi and 52- The systematic shift Ae^ s is an additional parameter 
which will be marginalized over. 

The ellipticity e of the total matter distribution is function of 
qi,q2, and of the orientation angles and if. The orientation angle 
6 e depends on the third Euler's angle, ip, too. The mass M200 and 
the concentration C200 of the total matter halo determine the lensing 
strength k s and the projected length-scale r a p. ft s and r s p depend 
also on the elongation, and, in turn, on qi, q2, 1? and if. 

Some a priori hypotheses on the cluster shape are needed to 
disentangle the intrinsic degeneracies. We applied some Bayesian 
methods already employed in either gravitational lensing investi- 
gations (|Oguri et al.|20"05||Corless et al.|2009 [ |Sereno et al.|2010| 
Sereno & Umets u|201 \\ or X-ray plus SZe analyses ( |Sereno et al.| 
2012 1 and extended them for our GL plus X-ray plus SZe analysis. 
We considered two kind of priors for the axial ratio of the mat- 



ter distributions. The distribution of gi obtained in high resolution 
iV-body simulations can be approximated as ( Jing & Suto 2002), 

p (gi )oce X p[- ^-f/ r "^ 2 l (14) 

where q^ — 0.54, <r s = 0.113 and 

r qi = (Af vir /Af*) ' 07nM(z) °"\ 

with M, the characteristic nonlinear mass at redshift 2 and A/ V i 
the virial mass. The conditional probability for 52 goes as 



(15) 



P(?l/?2|«l) 



2(1 



2<?l/<?2 



(16) 



for qi/q2 ?"min = max[qi, 0.5], whereas is null otherwise. 

We considered also a uniform distribution for the axial ratios 
in the range g m i n < qi ^ 1 and q\ ^92 1- Probabilities are 
defined such that the marginalized probability P(qi) and the condi- 
tional probability P(q2\qi) are constant. The probabilities can then 
be expressed as 



P(Qi) = V(l - <?m"0 

for the full range q m i n < q\ ^ 1 and 

pMgi) = (1- qiY 1 



(17) 



(18) 



for 52 ^ qi and zero otherwise. A flat distribution is also compati- 
ble with very triaxial clusters (qi < q2 <C 1), which are preferen- 
tially excluded by A-body simulations. We fixed g m i n = 0.1. 

For the axial ratio of the gas we considered a uniform distri- 
bution in the interval q± ^ q\ CM ^ 1- The prior on q\ CM is then 
similar to that of 92 in the case of the uniform distribution for the 
matter axial ratios. 

For the orientation, we considered a population of randomly 
oriented clusters with 



p(cosi}) — 1 

for ^ cos 1? ^ 1 and 

p(f) = - 



(19) 



(20) 



for -tt/2 sC f sC 71-/2. 

For the remaining parameters, we employed uniform distribu- 
tions. 



7 RESULTS 

Trends and correlations retrieved in a triaxial analysis are pre- 
dictable based on simple considerations. The efficiency of a halo 
as a lens is larger the larger either the mass or the concentration, 
all the rest being equal. The strength of a lens also increases with 
triaxiality and elongation along the line of sight, which boost con- 
vergence and cross section for lensing. A spherical analysis can 
not account for orientation and shape effects so that mass and con- 
centration estimates are systematically biased higher if the halo is 
oriented along the line of sight. On the other hand a full triaxial 
analysis is not affected by such bias, and the search for the best 
values of mass and concentration is not artificially constrained in a 
biased niche of the parameter space. 

The lensing analysis provides a reliable estimate of the pro- 
jected total mass. The projected ellipticity can be estimated either 
for the total mass from lensing or the ICM from X-ray maps. Meth- 
ods combining X-ray plus SZe constrain the elongation of the gas. 
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Priors 


M200 




(1 1 

y 1 






COS $ 


9i 


[lO 15 Af ] 












flat 


1.33±0.17 


7.8 ±0.7 


0.72±0.11 


0.8 


6±0.11 


0.57±0.29 


AT-body 


1.34±0.18 


7.6 ± 1.0 


0.51 ±0.07 


0.6 


9±0.12 


0.76 ±0.21 



Table 2. Intrinsic parameters of the matter distribution inferred assuming different priors on the axial ratios. Reported values are the mean and the variance of 
the posterior PDF. 



Flat q-distribution 




./V-body q-distribution 




Figure 2. PDFs for the matter halo parameters. Panels from the left to the right are for M200, C200, 91, 12 and cos respectively. The white binned histograms 
are the posterior PDFs. The dark grey histograms are the theoretical expectations given the measured mass in the left panel. For the concentration, we exploited 
the c(M) relation from |Duffy et al.H2008} . For the axial ratios, we considered predictions from ./V-body simulations ( Jing & Suto 20021. Intermediate grey 
bins are shared by the two plotted distributions. Thin lines in the q-panels denote a flat q-distribution. The thin line in the cos $ panel stands for random 
orientation. In the top row, we plot the inversion results under the prior hypotheses of flat q-distribution and random orientation angles. In the bottom row, 
prior assumptions are JV-body like q-distribution and random orientation angles. Masses are in units of 10 15 Mq. 



Mass, concentration and shape and orientation parameters have 
then to be chosen such that they reproduce the total lensing effi- 
ciency measured by GL within the geometrical limits which come 
from the knowledge of the size of the ICM along the line of sight 
and in the plane of the sky, as well as the size of the total matter in 
the plane of the sky. 

The resulting counterbalancing effects can be seen in Figs. [3] 
and|4] Lesser values of mass and concentration are compatible with 
halos elongated along the line of sight (cosi? < 1), i.e., elonga- 
tion supplies the lensing strength lost due to the decline in mass 
or concentration. Being the size in the plane of the sky fixed by 
observations, such elongation has to be fueled by a larger degree 
of triaxiality (lower values of 51). To account for the X-ray plus 
SZe constraints, the gas has to be more triaxial too (larger values of 
e ICM /e Mat ). 

The full triaxial analysis is required to quantify these trends 
and correlations. The results of our analysis are summarized in 
Table [2] They agree with our prev ious lensing {S ereno & Umetsu 
|201 1| | or X-ray plus SZe analysis ( |Sereno et al.|2012| . This is ex- 
pected since we used the same data sets and similar methodology. 
There are two main improvements: f) intrinsic halo parameters are 
derived with a better accuracy; ii) we rely on a considerably smaller 
number of hypotheses. In particular, the X-ray plus SZe data give 
directly the information on the elongation of the halo and we do not 
have to assume any orientation bias to reconcile results with theo- 
retical predictions. Likelihood functions combining X-ray plus SZe 
with lensing are more peaked than the corresponding functions ex- 
ploiting only one data-set and final results are less sensitive to the 



a priori hypotheses. Results assuming either a flat distribution of 
axial ratios or inputs from iV-body simulations are in very good 
agreement. Thanks to the information from X-ray plus SZe, which 
positively constrains the orientation, the alignment of the cluster is 
clearly seen whatever the priors. 

The posterior distributions were investigated by running four 
Markov chains and checking for convergence. The marginalized 
ID posterior probabilities are plotted in Fig. [2] whereas the bidi- 
mensional probabilities are represented in Figs. [3] and [4] Results 
for mass and concentration are sensitive to the assumed priors to a 
very small extent. The multi-probe approach is dominated by the 
data, and in turn by the likelihood, whereas the priors play a negli- 
gible role. 

As already touched in Sec. [3] our SL and WL constraints 
on the projected NFW parameters are marginally consistent and 
it makes sense to combine them. We used lensing data-sets al- 
ready exploited in previous wor ks (|Limousin et al.|2007| Umetsu| 
|& Broadhurst|[2008l |Coe et al.||2010| |Umetsu et al.||2011| > and 
we found consistent results. In particular, we refer to Sereno & 
|Umetsu| ( |2011| > for a detailed discussion of the compatibility of WL 
and SL results for A1689 in a triaxial context. For completeness, 
we quote here the results for our multi probe approach exploit- 
ing just one lensing data-set. Exploiting only the WL, we found 
M200 = (1.39 ± 0.32) x 10 15 M Q and c 20 o = 13.7 ± 3.9. SL 
alone favors less concentrated haloes, C200 
the mass is poorly constrained. 



4.8 ± 0.8, whereas 
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Figure 3. Contour plots of the bi-dimensional marginalised PDFs derived under the prior assumptions of uniform g-distribution and random orientation angles. 
Contours are plotted at fraction values exp(— 2.3/2), exp(— 6.17/2), and exp(— 11.8/2) of the maximum, which denote confidence limit regions of 1, 2 
and 3<r in a maximum likelihood investigation, respectively. In the M200 — C200 plane (top right panel), the thin and thick full line are the predictions from 
Duffy et al. 1 2008 1 or Prada et al. ( 20121, respectively. Dashed and long-dashed lines enclose the 1 and 3<r regions for the predicted conditional probability 
c(M) relation of |Duffy et al.||2008} , respectively. In the q\ — 92 plane, the thick full, long dashed and dashed lines limit the 1, 2 and 3cr confidence regions 
for the W-body like expected distributions. 



7.1 Mass and concentration 

The final results on mass and concentration are independent of the 
assumed priors. A 1689 is a very massive cluster with high concen- 
tration. Comparisons with theoretical predictions point to an over- 
concentrated cluster, see Figs.[2][3]and[4] We considered the c(M) 
relation from the analysis of a full sample of clusters in Duffy et al. 
(2008). The estimated concentration is in agreement with the tail 
at large values of the population of clusters of that given mass. Re- 
cently, Prada et al. 1 2012| > claimed that the c(M) relation features a 
flattening and upturn with increasing mass with substantially larger 



estimated concentrations for galaxy clusters. In that case, the agree- 
ment with theoretical predictions further improves. 

The distribution of the ellipsoidal radius r2oo can be derived 
too. We found r 2 oo = 3060 ± 290 kpc (2530 ± 280 kpc) for 
a priori TV-body (uniform) axial ratios. The spherical rf^, which 
comprises the same over-density and mass of the ellipsoidal T2oo 
and can be more useful for comparison with works based on spher- 
ical analyses, is 2135 ± 90 kpc in both cases. 

Concentration and orientation are strongly correlated, see 
Figs. [3] and [4] For TV-body axial ratios, the concentration corre- 
sponding to orientation angles 40deg < $ < 45deg is C200 = 
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Figure 4. Contour plots of the bi-dimensional marginalised PDFs derived under the prior assumptions of jV-body like axial ratios and random orientation 
angles. Contours and lines are as in Fig. [3] 



8.2±0.4.For0deg < & < 5deg, c 20 o = 6.2 ± 0.3. Even if the a 
posteriori probability distribution in the full parameter space peaks 
at lower concentrations and more pronounced alignments with the 
line of sight, the tail corresponding to larger values of the orienta- 
tion angle shifts the peak of the marginalized distribution towards 
larger concentration values. 



7.2 Shape and orientation of the matter distribution 

We found evidence for a triaxial matter distribution. This claim 
comes mainly from the observed ellipticity of the total projected 
matter. The measured value of e deviates from the spherical case, 
e = 0, by 2- and 5-a in the WL or SL analysis, respectively, see 



Table [T] The case of a prolate/oblate ellipsoid with the symmetry 
axis along the line of sight is ruled out at the same confidence level. 

When assuming a priori flat distribution, the posterior proba- 
bility of qi peaks at ~ 0.8 with a tail in correspondence of more 
triaxial shapes which brings the mean value at gi ~ 0.7. The large 
tail at small values makes the results fully consistent with theoret- 
ical predictions. When assuming the prior from iV-body simula- 
tions, the posterior probability of q 1 follows the prior, but with a 
marginal shift towards rounder values and a smaller dispersion. 

The intermediate to major axis ratio is less constrained. The 
prior plays an heavier role for the final distribution of q2, even if 
values of q-z ~ 0.7-0.9 perform well under different assumptions. 
Prolate configurations (q\ — qz) are slightly preferred over oblate 
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Priors 


q lCM 




e ICM y e Mat 


flat 


0.76±0.10 


0.87 ±0.10 


0.94±0.05 


JV-body 


0.56 ±0.08 


0.72 ±0.11 


0.96 ±0.04 



Table 3. Inferred intrinsic parameters for the shape of the gas distribution 
(axis ratios q\ CM and q\ CM ) and relation with the total matter distribution 
(e ICM /e Mat ). ql CM and e ICM /e Mat are derived parameters. Central val- 
ues and dispersions are the mean and the standard deviation of the PDF, 
respectively. 



mass one, see Table[T| even if not by a large margin. From the three- 
dimensional analysis, we found e ICM /e Mat ~ 0.95, see Table^ 
Very low values of e ICM / e Mat correspond to a nearly spherical gas 
distribution, which is excluded. On the other hand, very high values 
of e ICM / e Mat would mean that the gas follows the matter distribu- 
tion rather than the gravitational potential and are not excluded by 
data. 

We found q[ CM - qi = 0.05 ± 0.04 (= 0.04 ±0.04), for a N- 
body like (flat) prior. The probability that the difference qi CM — qi 
is larger than 0.1 is ~ 15 (6) per cent for a iV-body like (flat) prior. 



ones (<j2 = 1), but triaxial shapes give much better fits than axially 
symmetric haloes. 

The halo turns out to be elongated along the line of sight. Bi- 
ased orientations are favoured even if a priori the orientations were 
random. The a priori probability of a randomly oriented cluster to 
have < 45 deg is ~ 29 per cent. A posteriori such probability is 
~ 75 (42) per cent assuming iV-body like (flat) axial ratios. This 
result could be obtained only combining the lensing analysis with 
the information from X-ray and SZe. The pure lensing analysis in 
|Sereno & Umetsu ( 201 1 ) had to assume a priori a biased orienta- 
tion to get similar results on mass and concentration. Orientations 
in the plane of the sky or intermediate inclinations are compatible 
with less triaxial shapes and more massive and concentrated haloes, 
but their statistical significance is lower than for alignments along 
the line of sight. 



7.3 Shape of the gas distribution 

The X-ray plus SZe part of the inversion method directly constrains 
the size of the gas distribution along the line of sight and in the 
plane of the sky. The axial ratios for the ICM, q\ CM and ql° M , 
are then determined with a better accuracy than their counterparts 
for the matter distribution, q\ and 52- Results are summarised in 
Table|3] Posterior PDFs are plotted in Fig.[5] Final results are nearly 
independent of the priors. The estimate of the elongation e^ M for 
the gas and the exquisite accuracy in the measured ellipticity of 
the X-ray surface, e x , drive the final results on orientation. Even if 
there is some interplay with GL, ICM shape and cluster orientation 
are mainly determined by the X -ray plus SZe likelih ood. Results 



Sereno et al. 



(2012|Q The gas is 



are then very similar to those in 

mildly triaxial, gi° M ~ 0.6-0.8. A tail of the distribution extends 
to very low values, gi CM < 0.4, which are associated with very 
well aligned configurations (cost? < 1). 

Some trends between gas shape and matter halo parameters 
can be seen in Fig.Jel] For e ICM / e Mat < i i the total matter would 
be forced to follow the gas shape and orientations, which are well 
determined by the X-ray plus SZe part of the analysis. The matter 
would then be slightly rounder and very well elongated along the 
line of sight. As a consequence concentration and mass would be 
lower. 

We made the theoretical assumption that the gas is rounder 
than the matter distribution and data further supports this view. The 
ellipticity of the ICM, e x ~ 0.15, is lower than the projected total 



1 This analysis differs in the used priors from |Sereno et ah j2012| . In 
Sereno et al. 20121, a flat prior on q* CM was used. Here, q| CM ^ qi, 
so that a priori p(gj ) is not flat but peaks at q\ CM = 1. Being the anal- 
ysis dominated by the likelihood, this brings about only a small difference 
in the final results. 



7.4 Comparison with the spherical approach 

The deprojection method we employed was based on a minimum 
set of hypotheses. Since we did not assume hydrostatic equilibrium, 
the analyses of the ICM, based on X-ray and SZe observations, and 
that of the total mass, relying on lensing, are mostly independent. 
Their only tie is of geometrical nature, since we required the ICM 
to share the same orientation of the total mass, both in the space 
and in the plane of the sky. This can be seen in Eq. (12) , where the 
mass and concentration enter the likelihood only through the lens- 
ing part, whereas the X-ray and SZe contribution is only related to 
the ICM shape, and as a consequence, to the orientation of the mat- 
ter halo. Since the total lensing strength depends on the shape and 
orientation, the X-ray and SZe can then play a role in the overall 
properties of the halo. 

By assuming the spherical geometry, ellipticity and elongation 
of the total and gas mass distribution are fixed (e = 0, eA = 1). 
Therefore, the likelihood in Eq. jT3} is reduced just to the lensing 
part, with the ICM that can affect the mass reconstruction from 
lensing only through the constraints of the halo centroid, since the 
orientation is not anymore a parameter. Under these assumptions, 
M 200 = (1-26 ± 0.12) x 1O 15 M and c 20 o = 7.8 ± 0.2. As 
expected, the central values are compatible with the full triaxial 
analysis but have associated a statistical error that is small because 
it does not include the systematic part due to the relaxation of the 
assumption on the geometrical shape. 

Values of C200 ~ 6, which are fully compatible with the tri- 
axial analysis in the likely case of nearly alignment with the line 
of sight, are excluded assuming a spherical shape. As discussed in 
Sec. [7] the SL and WL analyses are marginally consistent assuming 
a triaxial form. The conflict between the results in the two lensing 
regimes is further aggravated in the spherical hypothesis. 



8 HYDROSTATIC EQUILIBRIUM 

A significant part of X-ray analyses relies on the assumption of 
hydrostatic equilibrium, 



VP-Tot = — PlCMV</!>Mat, 



(21) 



where Prot (= Prh±PnTh) is the total pressure, picm the gas den- 
sity and 0Mat the gravitational potential. When hydrostatic equilib- 
rium holds, the pressure is only thermal, Prh = fcBTViicM for an 
ideal gas, where fee is the Boltzmann constant. Neglecting a non 
thermal contribution P n Th, such as bulk and/or turbulent motions, 
systematically biases low the X-ray mass determination of the clus- 
ter (Meneghetti et al.|2010||Rasia et al.|2012| >. 

Assessing the level of hydrostatic equilibrium in a cluster 
can be problematic and usually relies on either multiple data sets 
jMorandi et al.| [bllbl or numerical simulations (Molnar et al. 
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Figure 5. Marginalized PDFs (plotted as white histograms) for the shape parameters of the gas distribution. Panels from the left to the right are for the 
axial ratios q\ CM , q\ CM and for e ICM /e Mat , respectively. The top (bottom) row refers to uniform (TV-body like) priors for the matter axial ratios. The 
grey histograms represent the expected distributions of the axial ratios for a gas in hydrostatic equilibrium under the derived mass density, i.e, by assuming 
e ICM /e Mat =0.7. 
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Figure 6. Marginalised PDFs in the plane of (/' versus matter halo parameters. Panels from the left to the right are for q\ -Mwo> -C200> -<?l> -Q2 
and -cosi9, respectively. Contours are as in Fig. [5] In the top row, we plot the inversion results under the prior hypotheses of flat (/-distribution and random 
orientation angles. In the bottom row, prior assumptions are TV-body like (/-distribution and random orientation angles. Masses are in units of 1O 15 M0. 



[20T0j >. Since we determined mass and shape of A1689 without re- 
lying on any assumption on the status of the cluster, we can check 
if and how A 1689 departs from hydrostatic equilibrium. Firstly, we 
know that gas in equilibrium follows the gravitational potential of 
the halo. In that ca se, gas and matter eccentricities are related and 

. This prediction can be com- 



„ICM / Jvlat 



0.7(Lee&Suto 



2003 



pared with our results, see Fig. |5] and Table [3] We found that the 
gas distribution is more triaxial than the shapes expected under the 
assumption of complete hydrostatic equilibrium. The distribution 
of e ICM / e Mat peaks at < 1, and significance at 0.7 is very low. 
The chance to have e ICM/ e Mat ^ q g is of order f < 1 per cent. 
We can conclude from this first test that the shape of gas and matter 
are only marginally compatible with the hypothesis of hydrostatic 
equilibrium, the gas being decisively more triaxial than expected. 



Radiative processes can make the gas more triaxial in the central 
regions ( |Lau et al.|201 1[ >, which would explain the high degree of 
gas triaxiality found in A 1689 

As a second test, we checked if the hydrostatic equilibrium 
condition expressed in Eq. {2TJ is fulfilled. To this aim, we recom- 
puted the posteriori probabilities for the cluster parameters under 
the sharp prior of e ICM / e Mat = 0.7. We exploited TV-body like 
priors for the axial ratios. The gravitational potential of the ellip- 
soidal NFW halo was computed using the approximated formulae 
in |Lee & Suto| ((2003 l. The final result is summarized in Fig. [7] 
which was obtained assuming that the gas follows the potential. 
Thermal pressure is systematically lower than what needed. Hy- 
drostatic equilibrium is compatible with our results at the 3a con- 
fidence level. The non thermal contribution to the total pressure 
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£ [kpc] 

Figure 7. Fraction of total contributed by the thermal pressure in A1689 
as a function of the ellipsoidal radius £ of the gravitational potential. The 
full thick line tracks the median; the dark (light) gray shadowed regions are 
limited by the 68.3 (99.7) per cent quantiles. The horizontal line at 1 is the 
value expected for an halo in hydrostatic equilibrium. Lines are dashed in 
the radial range where results are extrapolated. 



required for equilibrium, i.e., the value of Brot solving Eq. (21) , is 
of the order of 20-30 per cent in the outer regions at £ ~ 1.4 Mpc, 
and even higher towards the center (~ 20-50 per cent). 

High resolution cosmological simulations showed that there 
is a significant contribution from non-thermal pressure in the core 
region of relaxed clusters ( |Lau et al. 2009] |Molnar et~aT]|2010| >. 
|Molnar et al. (2010) found in ten simulated massive relaxed clus- 
ters that non-thermal pressure support from subsonic random gas 
motions can contribute up to 40 per cent in the inner regions and 
up to 20 per cent within one tenth of the virial radius. They also 
found that the non thermal contribution increases with radius in the 
very outer regions. |Lau et af7| (2009 1 found a similar level of non 
thermal pressure of the order of 5-15 per cent at about one tenth 
of the viral radius, also increasing with radius in the outer regions. 
These trends are retrieved in Fig.^] even if we have to caution that 
XMM spectroscopic data covers the cluster up to < 900 kpc so that 
the decrement at large radii is an extrapolation of fitted profiles. 



9 COMPARISON WITH PREVIOUS ANALYSES 

Our multi-probe analysis of A1689, combining a data-set spanning 
from X-ray to lensing to SZe data, can be compared to some re- 
cent works. Ett ori et aL] (2010] l performed an X-ray analysis and 
recovered the profiles of gas and dark mass in A 1689 under the 
assumption that hydrostatic equilibrium holds between the intra- 
cluster medium and the gravitational potential. Even if Ettori et al. 
(2010) used different techniques from ours, they still assumed a 
NFW functional form for the dark matter distribution and based 
their analysis on XMM data. Their results are then comparable to 
ours and any discrepancy should be interpreted in terms of either 
triaxiality or non-thermal pressure. They found C200 — 8.31^q g 3 , 
andM 20 o = (0.74 ±0.04) x 10 15 A/ Q for the dark matter compo- 
nent and caoo = 7.90i ;^ and M 20 o = (0.91 ± 0.05) x 10 15 M o 
for the total matter distribution. The estimated concentration is in 
excellent agreement with our findings whereas M200 is lower by 
-~30 per cent. This is exactly how much we expect the X-ray mass 
to be biased low due to the level of non-thermal pressure. Being 



A1689 nearly aligned with the line of sight, the spherical symme- 
try hypothesis brings about a bias high of ~5— 10 per cent ( |Gavazzi| 
|2005[ >, which is secondary with respect to the systematic error con- 
nected to the non-thermal pressure. 

Previous results from gravitational lensing found a concen- 
tration of C200 ~ 5-7 for strong lensing analyses or ~ 8-12 for 
weak lensing analyses. A collection from literature can be found 
in |Corless et al.| ( |2009| table 4) and |Coe et al.|(2010| table 2). 
Our results are in line with what already understood. In general, 
triaxiality and orientation issues reduces the gap between SL and 
WL analyses by allowing for smaller concentrations and larger dis- 
persions (Oguri et al.| 2005 , Corl ess et al.|20 09 , Sereno & Umetsu 



201 1). We retrieved the same trends already discussed in Sereno 



& Umetsu] (20 1 1| >, whose analysis is here potentiated in two main 
respects. First, we performed a non-parametric analysis of the SL 
core regions, which enables us to fully explore the parameter space 
and to be less sensitive to local minima. Secondly, the additional in- 
formation from X-ray plus SZe data reduces sensibly the parameter 
uncertainties and makes results much less sensitive to priors. 

|MorandietaT]pOTT) found C200 = 5.6 ± 0.4 and qx ~ 0.5 
from a combined lensing plus X-ray analysis. The level of triax- 
iality is similar to what found here whereas the concentration is 
quite smaller. This can be due to their hypothesis of halo perfectly 
aligned with the line of sight, which bias low the concentration. 

The level of hydrostatic equilibrium has been addressed in 
several works. [Morandi et al] (201 1| > found non-thermal pressure 
of order of ~ 20 per cent. They assumed the fraction to be con- 
stant with radius and fixed the gas shape to what expected in hydro- 
static equilibrium. |Lemze et al.| {2009 1 derived a temperature pro- 
file consistent with hydrostatic equilibrium from combined Chan- 
dra X-ray brightness measurements and joint strong/weak lensing 
measurements. They found that the resulting equilibrium tempera- 
ture exceeds the observed temperature by 30 per cent at all radii. 
The existence of significant cool and dense gas components was 
proposed as the source of this temperature discrepancy ( Ka wahara| 
|et al.|[2007) |Lemze et al.||2009) . Assuming spherical symmetry, 
|Peng et aL] (2009) derived an hydrostatic mass within the central 
200kpc/i _1 region of A1689 from Chandra observations, which 
is about 30-50 per cent lower than lensing-based mass estimates. 
They also showed that dense and cool gas clumps alone can not 
cause such a strong bias in the X-ray temperature determinations 
and proposed an orientation bias to reconcile mass determinations. 
Kawaharad a et al.|(2010) detected anisotropic temperature and en- 
tropy distributions in Suzaku images of the cluster outskirts corre- 
lated with large-scale structure of galaxies, with regions of low gas 
temperature and entropy deviating from hydrostatic equilibrium. 
They estimated that the thermal gas pressure within half the virial 
radius is at most about half of the equilibrium pressure required 
to balance the gravity predicted by gravitational lensing under the 
hypothesis of spherical symmetry. 

Our derived ratio of thermal to total pressure as a function of 
radius agree withMolnar et al. (2010), who assumed spherical sym- 
metry and used Chandra data for the temperature profile, which are 
larger by 10-20 per cent than the XMM measurements considered 
here and consequently may overestimate the thermal contribution. 
They found Pth/frot — 0.6 within the core region. 



10 CONCLUSIONS 

The multi-probe approach to galaxy clusters can tackle in an ef- 
ficient way one of the classical problems in astronomy, the deter- 
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mination of the mass. An unbiased estimate requires at the same 
time the knowledge of the shape and the orientation of the halo, 
of the concentration and of the equilibrium status of the gas. This 
complete picture can be achieved with an analysis exploiting gravi- 
tational lensing, which describes the total mass, and X-ray and SZe 
observations, which directly constrain the gas. 

We proposed a novel method based on minimal geometrical 
principles. We only assumed the gas and the total matter distribu- 
tions to be approximately ellipsoidal and co-aligned. The matter 
and the ICM were separately modelled with parametric profiles, 
suitable with comparison with numerical simulations and theoreti- 
cal predictions. We did not assume any hypothesis on the equilib- 
rium status of the gas or the profile of non-thermal pressure. In this 
regard, our method is alternative to other recently proposed multi- 
probe approaches (Mo randi et al.|20lT||01 lb| >. 

We obtained pixellated maps of the projected mass density 
covering a large radial range by combining strong lensing in the 
inner core and weak lensing in the outer regions. The maps were 
then fitted with a NFW profile. Photometric and spectroscopic X- 
ray data, as well as the SZ temperature decrement, were described 
with a single parameterization for the gas density and temperature 
profile. The combined X-ray plus SZe analysis enabled us to de- 
termine the elongation of the gas along the line of sight, the only 
geometrical quantity not defined in the plane of the sky that can 
be directly derived from the combined analysis of projected maps 
( |Sereno|2007| l and crucial information to constrain gas shape and 
orientation. All pieces of information were finally combined in a 
single Bayesian statistical analysis to infer the intrinsic parameters 
of the halo. 

The method was applied to A1689. We proposed the first 
multi-probe analysis of A1689 combining a data-set spanning from 
X-ray to lensing to SZe data. We could obtain an unbiased picture 
of the cluster. A1689 is massive, M 20 a = (0.9±0.1) x W 15 M e /h, 
and slightly over-concentrated C200 = 8 ± 1. The halo is triaxial 
(qi ~ 0.5±0.1) and aligned with the line of sight. The high degree 
of triaxiality (q\ CM ~0.6-0.8) of the gas distribution shows a de- 
viation from hydrostatic equilibrium, which would prefer rounder 
gas shapes. A significant contribution of non thermal pressure is 
required for equilibrium, ~ 20-50 per cent in the center and ~ 20- 
30 per cent in the outer regions. This level of non-thermal pressure 
support is consistent with what found by M olnar et al.| (2010) using 
a sample of massive relaxed clusters drawn from high resolution 
cosmological simulations and with recent findings in MS2137.3- 
2353 by |Chiu et aL| l |2012} , who found a 40-50 per cent contribu- 
tion in the core assuming a spherical model. 

The measurement of non-thermal pressure requires an unbi- 
ased knowledge of the cluster shape. [Chiu et al.| ( |2012| l found that 
the effect of the alignment of the major axis with the line of sight is 
to decrease the non-thermal pressure support required for equilib- 
rium at all radii without changing the distribution qualitatively. Dif- 
ferent counterbalancing factors can play a role in the determination 
of the cluster mass with X-ray methods, M%- Under the hypothesis 
of hydrostatic equilibrium, neglecting non-thermal processes, Mx 
is usually biased low by ~ 20-30 per cent (Rasia et al. 2012). Non 
thermal pressure would then make the lensing signal greater than 
expected given the X-ray derived mass. 

On the other hand, if a cluster aligned with the line of sight is 
considered spherical, Afx is biased high by < 5—10 per cent 
(Gavazzi 2005). The effect of the elongation on lensing is even 
more influential, since the central projected mass density of the lens 
is directly proportional to the extension of the halo along the line 
of sight. 



In this regard, our method seems particularly promising. 
Shape and inclination are measured and the orientation bias is over- 
come. The condition for hydrostatic equilibrium is not used to de- 
rive the mass and can be employed to determine the non thermal 
contribution to the pressure. Degeneracies are broken thanks to the 
joint multi-wavelength data sets giving a reliable picture of the clus- 
ter status and properties. 

The assessment of the over-concentration depend on the real 
presence of an upturn in the mass-concentration relation for high 
mass and redshift clusters (Prada et al. 2012). Massive systems are 
likely identified when they are substantially out of equilibrium and 
in a transient stage of high concentration (Ludlow et al. 2012). The 
upturn should disappear when only dynamically-relaxed systems 
are considered. Our results seem to support this picture. A1689 is 
massive and still not settled in hydrostatic equilibrium which pro- 
pounds an high value of concentration as observed. 
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